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Random percolation can be fully interpreted as a confining pure gauge theory. With numerical 
high-precision measurements of Polyakov-Polyakov correlators at finite temperature, we could 
well observe the presence of shape effects due to rough fluctuations of the confining string, in 
complete agreement with the universality predictions up to the next-to-leading order. 
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1. Introduction 

It is now widely believed that a D-dimensional confining gauge theory can be described, in 
the infrared limit, by a suitably chosen effective string theory. If the interquark separation is large, 
indeed, the relevant degrees of freedom are the oscillations of the confining string worldsheet in 
the D — 2 transverse directions to the loop plane. 

The Nambu-Goto action, which is proportional to the string worldsheet area, is however con- 
sistent under Lorentz covariance only in dimensionality D = 26; by slightly modifying the action 
with the addition of a nonpolynomial term, one gets the action ^ 

r 9,r. . „ n (dlX v d_X v )(d + X v d 2 X v ) - 
S = ajd 2 ^ {d + X^X» v {d+x ^ + xii y ) + 0(d- 3 )} , (1.1) 

which fixes the above problem and, while reproducing the Nambu-Goto features, gives physical 
predictions in any D>3. With the help of open-closed string duality in D = 3 [Q], and from some 
exact calculations available in any dimension, it has been argued that these effective theories are, 
to some finite order, universal. 

The validity of the effective string picture for confinement is strongly related to the phe- 
nomenon of roughening: it is exactly when the loop surface can undergo quantum fluctuations 
on any length scale that the identification with a bosonic, massless string is possible. 

By applying the Riemann zeta function regularisation prescription (to the Nambu-Goto effec- 
tive theory), it is possible to work out the functional form of observables such as Wilson loops in 
the form of an expansion in the small quantity 1 / (a - A), with a string tension and A the minimal 
area spanned by the loop contour [Slj. 

This work aims at determining the validity of the universality prediction by comparison with 
Monte Carlo data collected in a particular gauge theory, namely the random percolation model. In 
the following Sections, the observable of interest is presented in more detail, then a brief description 
of the gauge theory is given; afterwards, we describe the technique we used to obtain the data, and 
the results and conclusions that we drew. 



2. The Polyakov-Polyakov correlation function 

We focused on the behaviour of the Polyakov-Polyakov correlation function at finite tempera- 
ture in a (2 + 1)-D system. That is, the lattice is a L x x L y x L slice, with L x and L y large enough to 
represent the spatial extent and L = j the inverse temperature. We considered a couple of Polyakov 
loops orthogonal to the spatial direction and at a distance of R lattice spacings; the (connected) cor- 
relation function in this case is denoted with ( P(0)P*(R) ). 

Since the subleading terms in the correlator are a shape effect rather than exhibit a size de- 
pendence, we introduce the aspect ratio of the cylinder bounded by the loops and the sides of the 
lattice: 




(2.1) 
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The Polyakov-Polyakov correlator is then expected to follow the next-to-leading order (NLO) 
prediction 

< P(0)P*(R) >= , (2.2) 

where the functions v\ (Dedekind eta), E% and £4 (second and fourth Eisenstein functions) are 
defined by: 



n {T) = qvY[{\-cn, (2.3) 

n=\ 

OO 

E 2 ( Z ) = 1-24 £(71(71)9", (2.4) 

n=\ 

CO 

E 4 (z) = 1 + 240 £ 3 (n)q n ; (2.5) 



n=\ 



the functions C7j(n) here represent the sum of the i-tb powers of all divisors of n. This functional 
form for < P(0)P*(R) > is expected to be universal, to this order, for a wide variety of string 
functionals in the Nambu-Goto family [Q]. 

A direct numerical observation of these extremely fine corrections is, in most lattice gauge 
theories, unfathomably beyond any computational possibility. However, by choosing a particu- 
larly simple gauge model, a better numerical accuracy can be reached and the prediction can be 
adequately tested. This has been done on the three-dimensional Z2 gauge model in [Qj. 



3. The random percolation model 

In this work, the model we chose as laboratory is the random percolation model. By suitably 
defining the observables, indeed, it can be shown that this theory, although remaining somewhat 
peculiar, behaves exactly as expected in a ordinary confining pure gauge theory [g]. This interpre- 
tation is supported by a number of theoretical arguments (such as the center vortex picture for con- 
finement and the reformulation of the dual gauge system in terms of Fortuin-Kasteleyn clusters) as 
well as many numerical evidences concerning physical expectations (universal ratios, string tension 
scaling, glueball spectrum and so on). Nevertheless, the percolation model is somewhat particular, 
having a trivial partition function and gauge group (Z = 1 and G = {e}); indeed it can be thought 
of as the q — > 1 limit of the 9-state Potts model. 

The important aspect is that, as will be shown in the next Sections, even though the explicit 
formulation of the model does not involve strings at all, clear signals of a rough string behaviour 
are identified. It is already known that in this model a loop obeys the predicted functional form 
at least at the leading-order (LO); the question we address here is whether the model confirms the 
expected universality even at the NLO. 

In the (bond-)percolation model, each link of the lattice is independently set to on or off ac- 
cording to some fixed probability p, which plays the role of a coupling constant. Given a lattice 
geometry, there exist a critical value p c , corresponding to the sudden appearance of an infinite 
connected cluster; this is a second-order transition point and will represent the deconfinement uan- 
sition, which can be easily mapped to a finite-temperature transition at T = T c because the critical 
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threshold depends on the temperature (thought as the inverse extension of the lattice in the periodic 
direction). 

The starting point in this pure gauge framework is the definition of loop observables: given a 
loop with contour 7, the value of the associated Wilson loop (or any other object, such as a couple 
of Polyakov loops) is defined as: 

W(y) = 1 no cluster is topologically linked to the contour 7; (3.1) 

W(y)=0 <^ otherwise. (3.2) 

From this definition follows naturally an area/perimeter law for large enough loops, which distin- 
guishes between the confined (p > p c ) and deconfined (p < p c ) phases of the theory, providing in 
the former a well-defined string tension and /^-dependent critical temperature that scale according 
to 

a = S(p- Pc ) 2v , (3.3) 
T c = T (p- Pc ) v . (3.4) 

The value of a loop, moreover, is insensitive to changes in the configuration that do not alter 
the loop structure; this is a sort of gauge invariance of the theory, besides drastically reducing the 
computational effort required for numerical investigations. 

4. Methodology 

We worked on a cubic lattice of size 128 2 x L, where L is the periodic inverse temperature. 
For each choice of the occupation probability p, corresponding to some deconfinement temperature 
T c = 1 /L c , we took systems at various temperatures in the range y ~ T ~ T c . 

For each of these systems, we measured (P(0)P(R) } (the dagger can be dropped, since this 
theory deals only with real loop values), varying the distance between the two Polyakov lines from 
Ro to R max ; the data are unbiased by the spatial finiteness of the system up to at least half of the 
system size, and we safely chose Rq = 8,/? max = 50. 

The expectation value of a Polyakov loop couple is defined in terms of topological linking with 
the rectangle (periodic in one direction) which has the loops as boundary; once enough numerical 



accuracy is reached, then, it is possible to fit the data sets to Eq. [2.2| , extracting the string tension 
by looking at a plateau in the choice of the fit interval. 

As for the dependence of the string tension from the temperature T = 1/L, from the same 



Eq. 2.2 it is possible to find that, for asymptotically large R, 



here, the symbol cr(L) denotes the physical quantity that scales well with the temperature, while 
o is only a parameter in the fit that will be attempted to check this temperature scaling, a would 
represent the zero-temperature string tension (thought of as a function of p) if the above NLO 
formula were exact. 
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We chose to focus our attention on two values of P: at P\ = 0.272380 (critical for l/T = 6), 
we considered the temperatures in the range l/T = 7, . . . , 15, while at P%= 0.268459 (critical for 
l/T = 7), we examined l/T = 8, ... , 15. 1 

To reach an acceptable statistics, we collected data from 10 5 configurations for each value of 
P,L. For different blocks of 8 values of R we used independently-generated configurations. 

There is another signal that the theory is in the rough phase: the quantity 

m = «r>.. ts *zL (4.2, 

should be a universal ratio with no adjustable parameters, that is, it is supposed not to vary for 
different realisations of the transition (meaning different choices of p and the corresponding T c ). 

4.1 Algorithm 

Due to the particular nature of the random percolation model, each configuration can be gen- 
erated independently from scratch, by simply filling an empty lattice with links that are randomly 
switched on. The tricky part is the measurement of the topological linking of the resulting cloud 
with a given surface; to this end, the first thing to do is to "clean up" the configuration, getting 
rid of dead ends and simply-connecting bridges between loop structures. This is done once for the 
whole configuration. 

On this "minimal gauge" configuration, then, the loop is measured in all possible spatial po- 
sitions with the technique of reconstructing each time the clusters in the configuration (by means 
of the Hoshen-Kopelmann algorithm) keeping track of the crossings of the loop surface, to detect 
nonzero winding numbers. 



5. Data and results 



We tried to fit the measured correlation functions to Eq. 2.2: for temperatures far enough from 
the deconfinement point (L greater than 9 for P\ and 8 for P2), and large enough distances (R > 8 
and R > 9 respectively), the NLO formula works well and gives reliable plateau for the string 
tensions o (Fig. [j]). 



The quantities a thus obtained can be inserted in the asymptotic formula to reconstruct 
the physical quantity o{T) = o{L). However, since the expansion is truncated at (ff(L~ 4 ), the 
precision of the data allows to observe a residual dependence on L (Fig. ^). 

One can then look for the next term in the corrections, keeping in mind that the predicted uni- 
versality will not be valid any more beyond the NLO. So, we made, for this first model-dependent 
term, the Ansatz 

in which a new parameter, C, has appeared. 

By fitting the data sets with this functional form, we could well identify stable values both 

for the "true" zero-temperature string tension o(T = 0;P) and the coefficient C in the L~ 6 term 
(moreover, the two results for C are compatible, as one would have hoped): 



Note that the zero-temperature critical point is located at p c (0) — 0.248812 . 
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Figure 1: Zero-temperature string tensions a obtained with the NLO formula at P\ at different temperatures. 
The quantity R m \ n denotes the lower end of the fit interval, the upper one being fixed to 50. Note that, for 
large enough values of L, a plateau appears quite soon. 
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Figure 2: Temperature dependence, for Pi, of the physical string tension o(T) (squares) and of the fit 
parameter a (circles), from which the former is extracted. 



p 


a(T = 0) 


C 


Pi 


0.0126(1) 


296 ±5 


Pi 


0.00925(2) 


302 ±4 



As a circular check, we inserted this non-universal correction in Eq. [2^, and re-fitted the data 
sets: the plateaux for different temperatures (as long as they are not too close to the critical point) 
now coincide, confirming the estimates for o(P) (Fig. |3|). 



We could also confirm that the adimensional ratio f(t) of Eq. |4.2| , as a function of the reduced 
temperature, indeed does not depend of the choice of P, as expected in any confining theory in the 
rough phase (Fig. [|). 



6. Conclusions 



In this work the universality of the string behaviour up to the NLO has been numerically 
proven in a particular realisation of a confining gauge system in the rough phase. In the percolation 
model we have confirmed that the Polyakov-Ployakov correlation function follows the expected 
behaviour; moreover, we could identify the first model-dependent correction in the behaviour of 
a(L) with a stable coefficient 2 . 



"The coefficient C appears compatible with the integer 300, suggesting it comes, as expected, from some multiplicity 
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Figure 3: String tensions plateaux obtained with the 0{L 6 ) formula for (P(0)P(R) ), for the data at P\. 
Now the position of the plateau is temperature-independent. 
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Figure 4: Behaviour of the universal ratio /(f) for the two examined values of P. Note that the function 
does not fall to zero exactly at T = T c . 

The universality of /(f) confirms the validity of the string picture in this model (for large 
enough interquark separations). However, the function /(f) does not seem to drop to zero at exactly 
T = T c . This should be due to the fact that, approaching criticality, the algorithm used has to include 
all topological classes of configurations. A further analysis on this aspect could be carried on. 
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count. However, by lack of information about the small-scale behaviour of the theory, we could not delve deeper into 
this issue. 
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